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We extend the early time ordering theory of Cahn, Hilliard, and Cook (CHC) so that our gen- 
eralized theory applies to solid-to-solid transitions. Our theory involves spatial symmetry breaking 
(the initial phase contains a symmetry not present in the final phase). The predictions of our gener- 
alization differ from those of the CHC theory in two important ways: exponential growth does not 
begin immediately following the quench, and the objects that grow exponentially are not necessarily 
Fourier modes. Our theory is consistent with simulation results for the long-range antiferromagnetic 
Ising model. 



The early time dynamics of systems quenched into un- 
stable states is of considerable interest [H 0, S S B IE 
0, HI • The first effective theory to treat this process was 
developed by Cahn and Hilliard [l[ and by Cook [§] . The 
CHC theory applies to processes such as spinodal decom- 
position and continuous ordering Q and predicts that the 
early evolution of the equal time structure factor follow- 
ing the quench is characterized by exponentially growing 
Fourier modes. A primary assumption of the CHC the- 
ory is that the initial configuration following the quench 
is an unstable stationary point of the free energy . This 
assumption fails for solid-to-solid transitions. 

In this paper we introduce a generalized theory which 
describes the early time kinetics of spatial symmetry 
breaking transitions. We will show that the kinetics can 
be separated into two well defined stages for systems with 
effective long-range interactions. In the first stage sym- 
metry breaking fluctuations grow non-exponcntially. In 
the second stage the evolution crosses over to exponen- 
tial growth analogous to CHC. When the initial phase is 
a solid, we predict that the objects which grow exponen- 
tially are not Fourier modes. 

Binder || has shown that the CHC theory is valid 
only when the effective interaction range is large, R ^> 
1 llfJl. as has been confirmed in Ising model simulations 
[111 Il2j. There is evidence that many physical systems, 
such as polymers [f| and metals 13 1 have effective long- 
range interactions. It is therefore natural to develop our 
theory in the context of a long-range model. If the order 
parameter is expanded in powers of R~ d / 2 14|, we can 
separate the background (of order R =1) from the noise 
induced fluctuations (of order R~ d / 2 ). We will see that 
the CHC theory describes evolving fluctuations with a 
fixed background. 

Our early-time theory applies when the background is 
non-stationary. Such a background occurs for example in 
the early-time kinetics of a solid-to-solid transition where 
the lattice quickly contracts or expands before before sig- 
nificant symmetry change occurs. We will show that 
the evolution of the background is described by noise- 
less dynamics which explicitly preserves the rotational 



and translational symmetries of the initial phase. We 
say that the transition involves spatial symmetry break- 
ing if the initial phase contains symmetries not present in 
the final phase. When spatial symmetry breaking occurs, 
we will show that the background evolves to a stationary 
phase that is unstable with respect to symmetry break- 
ing fluctuations. Our theory predicts that the growth of 
these fluctuations changes from non-exponential (stage 1) 
to exponential (stage 2) when the background converges. 
The growth of the symmetry breaking fluctuations is de- 
scribed by a linear theory which is valid for a time of 
order lni? For R sufficiently large, both stages of 

our theory are observed. 

We develop our theory in the context of a time- 
dependent Ginzburg-Landau model with explicit long- 
range Kac interactions fl5| . The non-conserved field 
4>(x, t) plays the role of an order parameter and evolves 
according to the Langevin dynamics 



d<f>{S,t) = _ M _5F R [4>] 



dt 



<50(x, t) 



Mrj(x, t). 



(1) 



Fe\4>] is the free energy of the configuration <j) at time 
t, and R represents the effective interaction range. The 
Gaussian white noise fj(x, t) has zero mean and second 
moment (fj(x , t)fj(x' , t 1 )) = k B TS(t - t')8(x - x'). We set 
M = 1 corresponding to the rescaling of time, t — * t' = 
t/M. The drift term is given by 



' 5ct>{x) 



= / d d x' K R {x' )4>(x - x") + f{4>{x)) + h (2a) 



(2b) 



= (A fl * + /(0(f)) + fc, 



where Ar is a Kac potential of the form Ar(x) — 
R~ d A(x/R). The function / represents entropic forces, 
and h is an external field or chemical potential. By con- 
vention, /(0)|0 =o = 0. 

We scale all lengths by R so that Eq. (TT]) simplifies to 



du(r, t) 
dt 



SF[u] 
5u{r, t) 



(3) 
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where r = x/R, and 



u(r, t) = 4>(x, t) 
5F\u 



5u(f) 



(A*u)(f) + f(u(f)) + h. 



(4) 
(5) 



The parameter R in Eq. |T]) appears solely as a prefactor 
to the noise term. The term r)(r, t) = R d / 2 f]{x, t) repre- 
sents Gaussian white noise with zero mean and second 
moment (rj(r, t)i](f, t')) = k B T5(t - t')5(r - f"), which 
follows from the identity a~ d 5(x/a) = 5(x). 

The form of Eq. ([3]) suggests expanding u in the small 
parameter Rr d l 2 : 



= u (°) + R-WuW + R- d u^ + 



(6) 



We substitute Eq. ^ into Eq. ^ and obtain the dy- 
namical equations 

— — = i i = A*u' ' +/ u (0) ) + 7 



= CU {1) +T], 



where 



£V = A*V + /'(" (0) )^, 



(8) 



(9) 



and /'(it) = df/du. We remark that the nonlinear dy- 
namics of ?r ) in Eq. |(7J) is deterministic and decoupled 
from higher orders. The dynamics of is stochastic, 
linear, and depends on w^ ) through C. 

As we have mentioned, the CHC theory emerges as the 
evolution of u^' when is a stationary point of the 
free energy. Let us see how this works for a disorder-order 
transition occurring after a rapid quench from infinite to 
finite temperature and h = (recall that /(0) = 0). At 
t = the system is initially disordered so Eq. ((7]) has the 
trivial solution = for all time. With this solution, 
Eq. {§} can be solved in Fourier space, 

uW(M)=«W(M)e D *' 

+ / di'e^'-'VM'): (10) 



where D(jfe) = A(fc) + /'(w (0) = 0). The structure fac- 
tor S(k,t) = (\4>\ 2 ) /V can be calculated using Eq. (g}, 
thus reproducing the CHC theory. For spin systems the 
volume V equals the total number of spins because the 
lattice spacing is taken to be unity. 

We can easily determine the time scale for which the 
CHC theory is applicable. Equation ^ is meaningful 
when the neglected 0(R~ d ) terms are small. One re- 
quirement is that R~ d l' 1 vS 1 > <C ~ 1. The exponen- 
tial growth of it' 1 ' from Eq. (fTU|) suggests that the linear 
theory breaks down at a time t ~ In R 0, [12] • 

For many phase transitions (such as solid-to-solid) we 
need to consider the evolution of both and 



Equation dH) predicts exponential growth of u^\t) when- 
ever C is time independent, which from Eq. ([7]) occurs 
when SF/Su^°\x,t) — 0. In general, the initial con- 
dition u^°'(t = 0) will not be such a stationary point. 
We will show that, due to symmetry breaking, con- 
verges to an unstable stationary configuration u* . Cor- 
respondingly, C will converge to a time independent op- 
erator. This instability of u* means that C will have 
positive eigenvalues, corresponding to the unstable sym- 
metry breaking growth modes. 

Let Gi and G / represent the symmetry groups of ro- 
tations and translations of the initial and final phases 
respectively. We will now show that vS°\r, t) is invariant 
under Gi provided that the potential A(r) shares the ro- 
tational symmetries of Gi . Discretization of Eq. ([7|) with 
the time step At yields 



u 



(0) _ „(o) 

t+At ~ 



+ At(A*u{ 0) +.f(u) u> ) + h) 



(oh 



(11) 



If is invariant under Gi at time t, it is invariant at 
time t + At, which follows from the properties of the 
convolution operation (*) and the rotational symmetries 
of A. Induction establishes the Gi symmetry of for 
all t 0. 

How does evolve for a phase transition with sym- 
metry breaking? We see from Eq. ([7|) that F[it(°)] is 
non-increasing. Physically, F must be bounded from be- 
low, so we expect u' ' to converge to some configuration 
it*. This convergence occurs on a time scale independent 
of R. Symmetry considerations ensure that u* is not 
the stable phase: there exists a spatial transformation g 
which is in Gi but not in G/ (the symmetry breaking 
condition) under which the configuration u* and not the 
stable phase is invariant. Because u* is not the stable 
phase, we expect that u* is an unstable free energy sta- 
tionary point. Simultaneous to the evolution of u^°\ the 
dynamical noise induces symmetry breaking fluctuations 
in u' 1 ' which are of magnitude R~ d l 2 . These fluctua- 
tions are unstable, and if has converged, will grow 
exponentially for a time proportional to In R, analogous 
to the predictions of CHC. 

We conclude that spatial symmetry breaking phase 
transition kinetics can be decomposed into two stages: 

1. t < to'- Nonlinear evolution of toward u*, a config- 
uration of minimum free energy subject to symmetry 
constraints. The configuration u* is not the stable 
phase. The dynamical equation for u^ 1 ' is linear but 
has an explicit time dependence. Note that to is inde- 
pendent of R. 

2- to < t < mi?: To a good approximation has con- 
verged to u* . The linear theory of tS 1 ^ becomes anal- 
ogous to the CHC theory, and describes exponential 
growth of the unstable symmetry breaking modes. 

These two stages are illustrated in Fig. []Jb). In con- 
trast, there is no stage 1 process in the CHC theory, as 
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FIG. 1: (a) The CHC theory is applicable if the initial config- 
uration u(0) is a free energy stationary point. CHC describes 
the immediate exponential growth of Fourier modes, lasting 
a time t ~ In R. (b) For symmetry breaking transitions (e.g., 
solid-to-solid) the early time kinetics of u(t) has two stages. 
In the first stage the leading order contribution to u evolves 
deterministically and non-linearly to a symmetry-constrained 
(shaded plane) free energy minimum u* in a time t ~ 1. In the 
second stage, symmetry breaking modes grow exponentially 
for a time t ~ ln_R. (c) Without symmetry breaking (e.g., 
solid-to-fluid) the leading order contribution to u evolves de- 
terministically to the stable phase u(oo) in a time t ~ 1. 



10 5 



3 

O 
-r- 

CO 



10 J 
10 1 
10° 



* 



-rfixxxxxx? 



X 512 x 512 spins, R = 85 
+ 1024 x 1024 spins, R = 170 



0.5 



1 1.5 

Time (MCS) 



2.5 



FIG. 2: Evolution of the structure factor peak ((|0 max | 2 ) /V) 
for the fluid-solid transition in the antiferromagnetic Ising 
model following a critical (h — 0) quench. The CHC theory 
correctly predicts exponential growth, beginning immediately 
after the quench for this transition. 



illustrated in Fig. QJa) . 

Phase transition kinetics without spatial symmetry 
breaking, such as solid-to-fluid, are qualitatively differ- 
ent. Here it^ will evolve to u* but, unlike the symmetry 
breaking case, it* is the stable phase because the symme- 
tries of Gi are included in Gf (no symmetries are broken 
in the transition) . Note that all the interesting dynamics 
in this transition occurs through u^°\ which is indepen- 
dent of the noise. This process is illustrated in Fig. [IJc). 

Let us see how exponential growth arises in the second 
stage of a symmetry breaking transition by considering 
Eq. ([7]). Because £ is a real and symmetric linear oper- 
ator, it has a complete orthonormal cigenbasis and real 
eigenvalues 17f . The eigenvectors of C are Fourier modes 
only if is uniform. We can express the dynamics of 
vS 1 ' in the eigenbasis of C: 



du v 



(i) 



dt 



Cwvp + Vv = A„u« + Vv, (12) 



where v and A„ represent the corresponding eigenvectors 
and eigenvalues of C. The subscripts indicate eigenbasis 
components, for example u v = J d d rv( f)u( f) and C vv / — 
J d d fv(r)Cv' (f) = X v 8 vv i. The eigenvectors are normal- 
ized, and we can show that {rj v {t)ri v i{t')) = S vv 'S(t — t'). 

For times t > to the operator C is time independent 
and Eq. (fT2|) can be solved directly: 

«W(t) = 4 1) (io)e A " ( *- to) + f dfe^-^if). (13) 



The exponential growth of is apparent. We can ex- 
press it' 1 ' in the Fourier basis, 



(14) 



where v{k) is the Fourier representation of the eigenvec- 
tor v. If R is sufficiently large and there is a single largest 
eigenvalue \ v , then a single eigenvector v will grow ex- 
ponentially faster than all others. In this case, and at 
sufficiently large times, we can approximate 



( 1 \k,t)^v(k)\u£\t )e x ^~^ + dt'e x ^-^T, v (t') 



(15) 

We see that the exponential growth of the eigenvector v 
implies exponential growth of all the Fourier modes of 

(\u^ l \k, £)| 2 ^} oc e 2Xvt . These Fourier modes even- 
tually dominate all other contributions to the structure 
factor S = (\4>(k, t)\j / V, provided that the linear theory 

is valid (t < In R). 

We now compare our generalized theory to simulations 
of the 2D antiferromagnetic Ising model with a long- 
range square interaction. This model contains a dis- 
ordered fluid phase, as well as clump, and stripe solid 
phases 0,0. In the clump phase, localized regions of 
enhanced magnetization are arranged on a square lattice. 
In the stripe phase, regions of enhanced magnetization 
are arranged in periodic stripes. All fluid-to-solid phase 
transitions involve symmetry breaking, as do the tran- 
sitions between clump and stripe phases. In contrast, 
solid-to-fluid transitions do not involve symmetry break- 
ing because the uniform fluid phase contains all possible 
spatial symmetries. 

For the Ising model with > 1 a free energy func- 
tional F[(j)\ can be derived 18j. Rather than the dy- 



namics of Eq. we use single-spin flip Monte Carlo 
dynamics to simulate the system. At each update a spin 
is selected at random and flipped with the Glauber tran- 
sition probability p = (1 + e l3AE )~ 1 . Time is measured 
in units of Monte Carlo steps per spin (MCS). 
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FIG. 3: Evolution of the structure factor peak ((|0 m ax| 2 ) /V) 
for the fluid-solid transition following an off-critical (h ^ 0) 
quench. The initial growth has two stages, confirming the 
prediction of our generalized theory. The first stage is non- 
exponential and is independent of R. The second stage is 
exponential growth of the symmetry breaking modes (in this 
case Fourier modes) in analogy to the CHC theory. 



In Figs. [2] and [3] we display the peak of the struc- 
ture factor, S(k,t) — ^(k,t) 2 \ /V, for fluid-to-solid 

phase transitions following critical (h = 0) and off-critical 
(h = 0.8) quenches. In both cases the temperature is re- 
duced from T — oo to 0.05. The critical and off-critical 
transitions are described, respectively, by the CHC the- 
ory and our generalization. As predicted, the off-critical 
dynamics can be separated into two stages: initial non- 
exponential growth followed by an extended period of ex- 
ponential growth. The growth modes are Fourier modes 
for both types of quenches considered because the initial 
phase is disordered. 

In summary, we have shown that the CHC theory can 
be generalized to describe solid-to-solid transitions. The 
key ingredient of this generalization is spatial symme- 
try breaking. The predictions of our generalized theory 
differ from those of the CHC theory in two fundamen- 
tal ways: (1) the exponential growth of the symmetry 
breaking modes does not immediately follow the quench, 
and (2) these symmetry breaking modes are not gen- 
erally Fourier modes. We have performed simulations 
of the long-range antiferromagnetic Ising model for the 
off-critical fluid-to-solid transition, and have confirmed 
the existence of a transient stage preceding exponential 
growth of the structure factor. A future paper will show 
simulations confirming this theory in symmetry breaking 



transitions [l8j |. Finally, we point out that our theory 
does not apply in the presence of symmetry breaking de- 
fects in the initial conditions. 
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